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Abstract. While often the point sources in X-ray surveys are dominated by AGN, with the high 
sensitivity of modern X-ray telescopes such as Chandra and XMM-Newton normal/starburst galax- 
ies are also being detected in large numbers. We have made use of Bayesian statistics for both the 
selection of galaxies from deep X-ray surveys and in the analysis of the luminosity functions for 
galaxies. These techniques can be used to similarly select galaxies from wide-area X-ray surveys 
and to analyze their luminosity function. The prospects for detecting galaxies and AGN from a pro- 
posed "wide-deep" XMM-Newton survey and from future wide-area X-ray survey missions (such 
as WFXT and eRosita) are also discussed. 
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X-RAY SURVEYS 

Extragalactic X-ray surveys are a powerful tool to study important source popula- 
tions such as active galactic nuclei (AGN), clusters of galaxies, and, recently, normal 
galaxies JH]. In the bandpasses of Chandra and XMM-Newton, X-rays penetrate col- 
umn densities up to 10^-^"^^ cm~^, and are therefore efficient at detecting moderately- 
obscured "Compton-thin" AGN. Compton-thick AGN (AGN with column densities of 
> 10^^~^^ cm~^) can also be detected in X-rays when prominent scattered emission is 
present (typically of order of ^ 1% of the intrinsic emission) ^ and can be detected in 
hard (E > 10 keV) surveys [3]. Therefore X-ray surveys are essential for a complete cen- 
sus of AGN. Normal galaxies can now be detected in large numbers thanks to the high 
sensitivity of Chandra and XMM-Newton, however since they have low luminosities 
{Lx < 10^^ ergs s"^), they are more difficult to detect than AGN. 



X-rays from Normal/Starburst Galaxies 

It has been known since the early 1980s that the X-ray emission of normal and 
starburst galaxies (galaxies with very high star-formation rates) are correlated with both 
the star-formation rate and stellar mass of the galaxies [4]. The physics behind this is that 
high-mass stars evolve rapidly (on time scales of 10^~^ years), and in turn explode as 
supernovae (SN). Occasionally these SN are detected in X-rays, however more often 
they heat the ISM to X-ray emitting temperatures (i.e., T > 10^~^ K) and produce 
neutron stars and black holes in X-ray binaries. X-ray binaries where the companion 
is a high-mass star (high-mass X-ray binaries) also have short evolutionary time scales. 



Therefore hot ISM and high-mass X-ray binaries track the current star-formation rate. 
Low-mass X-ray binaries have longer evolutionary time scales (on the order of Hubble 
times), and therefore track the integrated star-formation history of galaxies (i.e., the total 
stellar mass). 



X-ray Galaxy Survey Strategies 

There are several approaches to surveying specific sources types with low fluxes such 
as X-ray observations of galaxies. Deep, pencil-beam surveys of course probe to the 
faintest fluxes however due to the limited survey volume tend to result in low numbers 
of rare, high-luminosity objects. Wide area surveys (e.g., XMM-COSMOS) detect large 
numbers of AGN but are usually too shallow or survey too small of an area to detect 
significant numbers of galaxies. Another approach is to correlate large catalogs (e.g., the 
RC3 catalog or the SDSS) with archival data. For example, ^400 galaxies have been 
detected in X-rays based on correlating the SDSS with the 2XMM catalog (see Parnau 
et al. these proceedings). We are pursuing this approach in the case of Chandra and 
XMM-Newton archival data by taking advantage of the XAssist pipeline processing of 
these datf] and correlating the fields with the RC3 catalog. Note that by working with 
the original data rather than simply a source catalog we will be properly integrating 
the X-ray flux over the full RC3 ellipse for each galaxy, which is often larger that the 
telescope PSF, and also will be able to compute upper-limits for galaxies not detected in 
X-rays. 

Finally we discuss the option of pursuing X-ray observations of statistically-complete 
samples that are not X-ray selected. We have been observing galaxies selected from 
the Nearby Field Galaxy Survey, which have very-well determined star formation rates 
[El]. We have received six XMM-Newton datasets and also included one serendipitous 
Chandra observation [6]. The X-ray/SFR correlation based on correlating the Chandra 
archive with the Kauffman SDSS galaxy catalog [7] is shown in Figure [HS], with the 
X-ray NFGS data also plotted. The NFGS points are consistent with either the lower 
X-ray/SFR normalization implied by the X-ray detected SDSS galaxies, or a break in 
the X-ray/SFR correlation. Clearly a larger unbiased sample is needed. 

BAYESIAN SELECTION AND ANALYSIS OF NORMAL 
GALAXIES IN DEEP X-RAY SURVEYS 

In §\ and M\, normal/starburst galaxies were selected from the Chandra Deep Fields 
(CDF) using a Bayesian model selection methodology. Our motivation for employing 
this technique was to directly take into account the X-ray measurement errors which can 
be significant for X-ray sources at the flux limit of the survey. 




FIGURE 1. The X-ray/SFR correlation based on correlating SDSS galaxies and the Chandra archive 
JH], with data from XMM observations of the Nearby Field Galaxy Survey (NFGS) added. 



Galaxy Selection 

The most discriminating features for separating galaxies from AGN were found to be 
X-ray hardness ratio, (H-S)/(H+S), where H and S are the numbers of photons above and 
below 2 keV, respectively, X-ray luminosity, and X-ray/optical flux ratio. In our analysis 
of GOODS data we also included the X-ray/near-IR flux ratio. We determined "parent" 
distributions for these parameters by selecting a sample of normal galaxies, type-1 AGN 
and type-2 AGN from the CDF South based on high quality optical spectroscopy, and 
then taking the mean and standard deviation (cj) for each parameter. 

The posterior probability for observing the parameters 0, where here 6 = 
{H R ^log Lx ^log Fx / Fopi^log Fx / F^i^} given the data D is given by Bayes' theorem: 
Pm{0\D) = pm{0)pm{d\6) / pm{D). The M subscripts denote that this is assuming a 
given model M, here galaxies, type-1 AGN (AGNl) and type-2 AGN (AGN2). If mul- 
tiple models are being considered, then the prior probability for each model must also 
be included. pmO are the "prior" distributions for the parameters for a given model M, 
for which we used the parent distributions discussed above. pm{D\G) is the likelihood 
function of observing the parameters 6 given the data. Often pm{D) is considered to 
be a normalization constant since it does not depend on the model parameters 0, and 
is defined to be pm{D) = / d6pM{0\D). However pm{D) is also often (perhaps more 
precisely) considered to be the marginal likelihood or Bayesian evidence for the model 
M. The relative probability of two competing models given the data is then given by the 



"Bayes Factor" or 

PM,{D) ^ !depM,{e\D) 

PM,{D) JdOpMMD) 

Sources with /^galaxy/ Pagni > 1 ^nd /^galaxy/ Pagni > 1 wcrc selected as galaxie^. Here 
we are assuming a flat prior on the numbers of sources for each model M, in other words 
we are assuming that the "true" number of normal/starburst galaxies, type-1 AGN and 
type-2 AGN are approximately equal. 



Evolution 

Evolution was observed qualitatively between the redshifts of 0.25 and 0.75 consistent 
with pure luminosity evolution, L* (z) = L* (z = 0) ( 1 + zY with ^ 3 by comparing the 
X-ray luminosity functions (XLFs) with the far-infrared luminosity functions [9]. Sub- 
sequently using GOODS data we fit the X-ray luminosity functions by using Markov- 
Chain Monte Carlo (MCMC) [10]. MCMC analysis results in a distribution of parameter 
values (the "chains"). This allows the direct visualization of posterior probabilities for 
important quantities, such as the change in L* between the low and high redshift XLFs 
(Figure[2l). Another key advantage of the Bayesian approach is that "derived" quantities 
such as luminosity density can be computed directly from the chain parameter values, 
allowing the posterior probabilities for these quantities to be visualized or summarized 
(i.e., with the mode as a "best-fit" value and the 68% confidence interval as the "er- 
ror bar") without questionable propagation of error MlUD. In future work we will be 
incorporating the additional 1 Ms of CDF-S data, as well as improving our MCMC anal- 
ysis of the XLFs to also include VLA radio and Spitzer mid-IR data (both of which are 
star-formation rate indicators and will help discriminate galaxies from AGN). 



WIDE-AREA X-RAY SURVEYS 

XMM Wide-Area Surveys 

While both XMM-Newton and Chandra are very sensitive X-ray telescopes, XMM- 
Newton has a larger field-of-view (FOV) than Chandra and is therefore more adept at 
wide-area surveys. While the total solid angle of 2XMM is large (hundreds of square de- 
grees), its coverage is non-uniform and potentially biased since the field selection is not 
random. XMM-COSMOS is a wide-area survey covering 2 deg.^ at an exposure of 40 ks 
per field, or a limiting point-source sensitivity of ^ 5 — 10 x 10~^^ ergs s~^ cm~^[12]. 
Here we briefly discuss a proposal to extend XMM-COSMOS out to 10 deg.^ at a sim- 
ilar mean exposure per fielco This level of sensitivity would be sufficient to detect the 
(scattered) X-ray flux from Compton-thick AGN with a spectral energy distribution sim- 



^ In (Td], a more conservative sample was selected with Bayes factors > 3. 
^ Submitted in the XMM AO- 8 proposal round, PI David Alexander 
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FIGURE 2. Posterior probability for the change in logL* between the z ~ 0.25 and z ~ 0.75 for early- 
type (left) and late-type (right) galaxies. The solid (black) lines show the posterior probabilities, while 
the dashed (red) lines show Gaussian distributions with the same mean and standard deviation as the 
posteriors, and the dotted (green) lines show the prior distributions (~ flat in these cases) 

ilar to NGC 6240 at z - 0.5 - 2.0 (see Figure[3l). This survey should result in 6000-8000 
AGN being detected, ^ 1000 with at least 100 photons (sufficient for crude X-ray color 
analysis) and ^ 400 with at least 300 counts (sufficient for spectral analysis, including 
the detection of Fe-K lines). This survey would also detect ^ 150 — 200 AGN at z>3 
15 — 20 at z>4). We expect up ^ 300 normal/starburst galaxies would be detected, 
based on the Ranalli et al. logN-logS[13]. The field selection is in the Spitzer-SWIRE 
area, and the proposed X-ray data along with large amount of ancillary data from other 
wavebands available in these fields will allow us to study the coeval evolution of AGN 
and star formation over a wide range of redshift, environment and luminosity. 



Future Missions 

There are several proposed future missions dedicated to X-ray surveys. eRosita (ex- 
tended Roentgen Survey with an Imaging Telescope Array) is an approved mission ex- 
pected to fly aboard the Spectrum X-Gamma Mission, although the launch date appears 
to be uncertaiii^. eRosita will perform several surveys, including a nearly all-sky shallow 
survey. The Wide Field X-ray Telescope (WFXT) is a proposed medium-class NASA 
mission to similarly perform several surveys. The limiting flux and survey solid angle 
for both the WFXT and eRosita surveys are plotted in Figure |4l The numbers of AGN 
and clusters expected to be detect by eRosita and WFXT are shown in Figure [5l We ex- 
pect on the order of '--^ 10^ normal/starburst galaxies to be detected in the eRosita surveys 
while ^ 10^ galaxies should be detected in the WFXT surveys. Clearly either mission 



http://www.mpe.mpg.de/proj ects .html#erosita 
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FIGURE 3. The SED of NGC 6240 (dashed curves, for redshifts of 0.5 and 2.0 and assuming a 
Lx = 5 X 10"^^ ergs s~^) shown with the Hmiting fluxes of the SWIRE survey (blue), all-sky sureys (red) 
and the proposed XMM-Newton survey marked. 
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FIGURE 4. The limiting flux and solid angle of the surveys from the future missions WFXT and 
eRosita, where both the sensitivity to AGN (point sources) and clusters are shown, along with other X-ray 
surveys. 



would drastically increase the numbers of X-ray detected sources and would be revo- 
lutionary. The high numbers of sources expected from WFXT is due to both a higher 
effective area and a smaller PSF (half-energy width of ^ 7" for WFXT compared to the 
field-averaged PSF of ^ 25 — 30" in the case of eRosita). WFXT also have the advantage 
of using a wide-field optical design [14], giving more uniform PSF and response across 
the 1 deg. FOV. 




FIGURE 5. Expected numbers of AGN (left) and cluster of galaxies (right) to be detected by WFXT 
and eRosita. Also shown are the number of clusters of galaxies expected to be detected by the South Pole 
Telescope (SPT) via the Sunyaev-Zel'dovich effect. 
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